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ABSTRACT 

Using data from the Spitzer Space Telescope Legacy Science Program "For- 
mation and Evolution of Planetary Systems" (FEPS), we have searched for debris 
disks around 9 FGK stars (2-10 Gyr), known from radial velocity (RV) studies to 
have one or more massive planets. Only one of the sources, HD 38529, has excess 
emission above the stellar photosphere; at 70 fim the signal-to-noise ratio in the 
excess is 4.7 while at A < 30 /im there is no evidence of excess. The remaining 
sources show no excesses at any Spitzer wavelengths. Applying survival tests to 
the FE PS sample and the results for the FGK survey published in Bryden et al. 



(120061 ). we do not find a significant correlation between the frequency and prop- 
erties of debris disks and the presence of close-in planets. We discuss possible 
reasons for the lack of a correlation. 
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Introduction 



During the last two decades, space-based infrared observations, first with IRAS and then 
with ISO and Spitzer, have shown that many main sequence stars are surrounded by dust 
disks (a.k.a. debris disks). These disks are generally observed by their infrared emission in 
excess over the stellar photosphere, but in some cases the disks have been spatially resolved 
and extend to 100s of AU from the central star. Dust particles are affected by radiation 
pressure, Poynting-Robertson and stellar wind drag, mutual collisions and collisions with 
interstellar grains, and all these processes contribute to make the lifetime of the dust particles 
significantly shorter than the age of the star/disk system. It is therefore thought that 
this dust is being re plenish ed by a reservoir of undetected dust-producing planetesimals 
(Backman & Paresce Il993l ). like the asteroids, the Kuiper Belt Objects (KBOs) and the 
comets in our solar system. This represented a major leap in the search for other planetary 
systems: by 1983, a decade before extra-solar planets were discovered, IRAS o bserva tions 
proved that there is planetary material surrounding nearby stars (Aumann et al. Il984l ). 



Preliminary results from Spitzer observations of FGK (solar-type) stars indicate that 
the frequency of 24 /zm excesses (tracing warm dust at asteroid belt-like distances) decreases 
from ~30%-40% for ag es <50 Myr, to ~9% for 100 Mvr- 200 M yr, and ~1.2% f or ag es >1 



Gyr (Siegler et al. 12006 



Gorlova et al. 2006; Stauffer et al. 2005: Beichman et al. 2005; Kim 



et al. 120051 ; Bryden et al. 120061 ). On the other hand, Bryden et al. (120061 ) estimated that 
the excess rate at 70 /iin (tracing colder dust at Kuiper Belt-like distances) is 13±5% and 
is not correlated with stellar age on Gyr timescales. It is also found that FGK stars show 
large variations in the amount of excess emission at a given stellar age, and that the upper 
envelope of the ratio of the excess e missio n over the stellar photosphere at 24 /im decays as 
l/t for ages >20 Myr (Siegler et al. l2006h . 



These observations are consistent with numerical simulations of the evolution of dust 
gen erated from the collision of planetesimals around solar-type stars (Kenyon & Brom- 
ley l2005j). These models predict that after 1 Myr there is a steady l/t decline of the 24 
jum excess emission, as the dust-producing planetesimals get depleted. It is also found that 
this decay is punctuated by large spikes produced by individual collisional events between 
planetesimals 100-1000 km in size. These events initiate a collisional cascade leading to 
short-term increases in the density of small grains, which can increase the brightness density 
of the disk by an order of magnitude, in bro ad agr eement with t he hig h degree of debris disk 
variability observed by Spitzer (Rieke et al. 120051 ; Siegler et al. 120061 ). 



However, these models do not include the presence of massive planets, and the study 
of the evolution of the Solar System indicates that they may strongly affect the evolution 
of debris disk. There has been one major event in the early Solar System evolution that 
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likely produced large quantities of dust. Between 4.5 Gyr and 3.85 Gyr ago there was a 
heavy cratering phase that resurfaced the Moon and the terrestrial planets, creating the 
lunar basins and leaving numerous impact craters on the Moon, Mercury and Mars. This 
"Heavy Bombardment" ended abruptly ~3.85 Gyr ago, and since then the impact flux 
has been at least an order of magnitude smaller. During the last 20-200 Myr of the "Heavy 
Bombardment" epoch, a period known as the "Late Heavy Bombardment" (LHB), there was 
an increased cratering activity that came after a relatively calm period of several hundred 
million years, that could have been created by a sudden injection of impact objects from 
the main asteroid belt into the terrestrial zone. The orbits of these objects became unstable 
likely due to the the orbital migration of the giant planets which caused a resonance sweeping 
of th e aster oid belt and a larg e scale ejection of asteroids into planet-crossing orbits (Strom 
et al. 120051 ; Gomes et al. 120051 ). This event, triggered by the migration of the giant planets, 
would have been accompanied by a high rate of asteroid collisions, and the corresponding 
high rate of dust production would have caused a large spike in the warm dust luminosity of 
the Solar System. Similarly, a massive clearing of planetesimals is thought to have occurred 
in the Kuiper Belt. This is inferred from estim ates o f the total mass in the KB region , 
30-55 AU, ranging from 0.02 M e (Bernstein et al. l2004h to ~0.08 M ffi (Luu & Jewitt I2OO2J) , 
insufficient to have been able to form the KBOs within the age of the Solar System (Stern 
1996). It is estimated that the primordial KB had a mass of 30-50 M ffi between 30-55 AU, 
and was heavily depl eted after Neptune forme d and started to migrate outward (Malhotra, 
Duncan and Levison I200CH : Levison et al. 120061 ). This resulted in the clearing of KBOs with 
perihelion distances near or inside the present orbit of Neptune, and in the excitation of the 
KBOs' orbits. The latter increased the relative velocities of KBOs from 10 m/s to >lkm/s, 
making their collisions violent enough to result in a significant mass of the KBOs ground 
down to dust and blown away by radiation pressure. 

The evolution of debris disks may therefore be strongly affected by the presence of 
planets: it its early history, a star with planetary companions may be surrounded by a 
massive debris disk while the planets are undergoing orbital migration, whereas at a later 
stage the star would harbor a sp arse du st disk after the dynamical rearrangement of the 
planets is complete (Meyer et al. l2006al ); at very late stages, 2-10 Gyr, the production of 
dust may undergo occasional bursts due to major collisions of planetesimals stirred up by 
the planets. In addition to their effect on the dust production rates, massive planets can also 
affect the dynamics of the dust grains. Examples include the trapping of dust parti cles in 



mean motion resonances a nd th e ir ejec tion due to gravitational scattering (Liou & Zook ll999 
More-Martin & Malhotra [200% I2OO5J) . 



In this paper, we search for debris disks around nine stars known from radial velocity 
studies to harbor one or more massive planets. These stars are drawn from the Spitzer 
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Legacy program "Formation and Evolution of Planetary Systems" (FEPS). The properties 
of the stars and their planetary companions can be found in Tables 1 and 2. The observations 
and data reduction are briefly described in §2 and the resulting spectral energy distributions 
(SEDs) are presented in §3. In §4 we explore the correlation of the frequency of dust emission 
with the presence of known planets by applyi ng su rvival tests to the FEPS sample and the 
FGK-star survey published in Bryden et al. (120061 ). Finally, §5 discusses the interpretation 
of our results. HD38529, the only planet star in the FEPS sample with an excess emission, 
is discussed in detail in Moro-Martm et al. (in preparation). 



2. Observations and Data Reduction 



An overview of the FEPS program is given in Meyer et al. (120041 . 12006b ). and a detailed 



description of the data acquisition and data reduction in Hines et al. (120051 ) and Carpenter 



et al. (in preparation). The Multiband Imaging Photometer for Spitzer (MIPS; Rieke 



et al. 120041 ) was used to obtain observations at 24 /zm and 70 jira using the small field 
photometry mode with 2-10 cycles of 3 and 10 s integration times, respectively. The data 
was first processed by the Spitzer Science Center (SSC) pipeline version S13, and further 
processing was done by the FEPS team, the details of which can be found in Carpenter et 
al. (in preparation). 

At 24 jum, point- spread-function (PS F) fitt ing photometry was performed using the 



APEX module in MOPEX (Makovoz et al. 120051 ) using a fitting radius of 21 pixels on the 



individual Basic Calibrated Data (BCD) images. Fluxes were computed by integrating the 
PSF to a radius of 3 pixels, and then applying an aperture correction of 1.600 to place the 
photometry on the same scale described in the MIPS Data Handbook. The S13 images were 
processed using a calibration factor of 0.0447 MJy sr" 1 . We adopt a calibration uncertainty 
of 4% as stated on the SSC MIPS web pages. 

The raw MIPS 70/im images were processed with the SSC pipeline version S13. The 
individual BCD images were formed into mosaics with 4" pixel sizes using the Germanium 
Reprocessing Tools (GeRT) software package S14.0 version 1.1 developed by the SSC. The 
GeRT package performs column filtering on the BCD images to remove streaks in the BCD 
images, and then performs a time median filter to remove residual pixel response variations. 
A 40" x 40" region centered on the source position was masked when computing the time 
and column filtering such that the filtering process is not biased by the source. The filtered 
images were formed into mosaics using MOPEX. Aperture photometry was performed on the 
MIPS 70/im mosaics using a custom modified version of IDLPHOT. The adopted aperture 
radius of 16" was chosen to optimize the signal to noise for faint sources. The sky-level was 
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computed as the mean value of the pixels in a sky-annulus that extends from 40" to 60". The 
photometry uncertainty is given by a = ^lcr sky y^7i^r] sky r] corr ^/l.O + N ap /N sky , where a sky is 
the standard deviation in the sky annulus surface brightness, Q is the pixel solid angle, N sky 
and N ap are the number of pixels in the sky annulus and in the aperture, and rj sky and rj corr 
are correction factors that account for the presence on the mosaic of non-uniform noise and 
of correlated noise, respectively. We used rj sky = 2.5 and r\ COTr = 1.40 (see full description in 
Carperter et al. in preparation). The adopted calibration factor is 702 MJy sr _1 / (DN s _1 ) 
with an uncertainty of 7% as described on the SSC MIPS web pages. 



The Infrared Spectrograph (IRS; Houck et al. 120041 ) was used to obtain low- resolution 
(R = 70-120) spectra from 7.4 /iin to 38 /mi, with integration times per exposure of 6 s and 
14 s for the Short-Low (7.4-14.5 /mi) and Long-Low (14.0-38.0 /mi), respectively. The data 
was initially processed with the SSC pipeline S10.5.0, with further processing described in 
Bouwnman et al. (in preparation). From the spectra, synthetic photometric points were 
calculated at 13 /im with a rectangular bandpass between 12.4/im and 14.0/im, at 24/mi 
with the same bandpass shape as the MIPS 24 filter, and at 33 /mi with a rectangular 
bandpass between 30 /mi and 35 /mi. The estimated calibration uncertainty in the synthetic 
photometric is 6% (Carpenter et al. in preparation). The spectra are generally not reliable 
beyond 35 /mi, although we found that for HD 6434, HD 121504 and HD 80606 it is very 
noisy beyond 34 /mi, 33 /mi and 30/iin, respectively, making the 33 /iin photometric points 
unreliable for the last two sources. The latter is flaged by the SSC as non-nominal, possibly 
due to a failure in the peak-up. 



The Infrared Array Camera (IRAC; Fazio et al. 120041 ) was used to obtain observations 
at 3.6 /im, 4.5 /im and 8.0 /im in subarray mode. Initial processing of the data was done 
with the SSC pipeline S13, with further processing as described in Carpenter et al. (in 
preparation). Aperture photometry on individual IRAC frames were performed using a 
custom modified version of IDLPHOT using an aperture radius of 3 pixels (1 pixel ~ 1.2"), 
with the background annulus extending from 10 to 20 pixels centered to the star. The 
internal uncertainty was estimated as the standard deviation of the mean of the photometry 
measured at the four dither positions. We adopted calibration factors of 0.1088, 0.1388, 
and 0.2021 MJy/sr per DN/s fo r IRA C 3.6, 4.5/mi and 8/mi respectively and calibration 



uncertainties of 2% (Reach et al. 120051 ) 



3. Spectral Energy Distributions and Excess Emission 



The SEDs are shown in Fig. 1 and include the Spitzer photometric measurements, 
observations made by IRAS, Tycho and 2MASS and, in some cases, upper limits at 1.2 
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mm from Carpenter et al. (120051 ). For each star, the results from the Spitzer photometric 
measurements and their internal uncertainties are listed in the first entry of Table 3 (in rows 
indicated by "obs"). The reported fluxes arise from both the photosphere of the star and 
the thermal emission of the dust (if present). For all targets, observations are sufficient to 
detect the photosphere of the star at all Spitzer wavelengths <33 /zhi, making it possible to 
detect small dust excesses (limited mainly by the calibration uncertainties). To estimate the 
contribution from the dust alone, we need to subtract the photospheric emission, given in 
the second entry of Table 3 (under "model"). The Kurucz model calculations are described 
in Carpenter et al. (in preparation). Infrared excesses can also be identified empirically from 
the color-color diagrams in Fig. 2, showing a narrow distribution of the ratio of the fluxes 
at 24 (Ltm and 8 /im (F24/Fg) compared to a wide distribution of the ratio of the fluxes at 70 
jum and 24 /mi (F 70 /F 2 4). This indicates that the flux at 24 /im is mainly photospheric and 
that the best indicator of the presence of a debris disks is the 70 /im excess emission. 

HD 38529 has the only robust detection of an excess at 70 /im with a signal-to-noise 
(SNR) in the excess of 4.7, a small excess at 33 /im (with a SNR in the excess of 0.7) and no 
excess <30 /im. Because of the slope of the spectrum across the MIPS 70 /im band, the 70 /zhi 
flux needs to be color-corrected by dividing the observed flux by 0.893 (assuming the emission 
arises from cold dust emitting like a blackbody at 50 K - see MIPS Data Handbook). This 
increases the 70 /im flux from 75.3 mJy to 84.3 mJy The centroid positions of the object in 
the 24 /im and 70 /im images are RA=05:46:34.88, DEC=+01: 10:04.61 and RA=05:46:34.79, 
DEC=+01:10:04.62, respectively, in agreement with the 2004.7 2MASS coordinates for HD 
38529 (RA=05:46:34.895 and DEC=+01:10:04.65; accounting for the proper motion of the 
star), and with the absolute pointing knowledge: better than 1.4" and 1.7" (1-er radial) at 
24 /im and 70 /im, respectively (Spitzer Observers Manual). Inspection of the images show 
that the 24 /im and 70 /im source is free of nearby point sources, and there is very little 
structure from galactic cirrus. Finally, it is unlikely that the emission at 70 /im comes from 
a bac kgrou nd galaxy within 2"of the st ellar p osition: from the background counts in Dole 



et al. l2004j and following Downes et al. (119861 ) we estimate a probability of 1.5x10 5 for 50 



mJy and 7.4xl0 -6 for 100 mJy. We therefore conclude that the observed 70 /im emission 
comes from HD 38529. Even though it is difficult to identify statistical trends from one 
detection, it is interesting to note that HD 38529 is the most luminous, most massive, and 
most evolved of the planet bearing stars in Table 1. Assuming V = 5.95 (Johnson 1966), a 
Hipparcos distance of 42 pc and no reddening, the object has an absolute visual magnitude 
of M„ = 2.81 and Log(L/L ) = 0.82, putting the star on the Hertzsprung gap, so it is clearly 
post-main sequence. 



In summary, we find that only 1 out of 9 of the planet-bearing stars show evidence of a 
debris disk. In the next section, we explore whether or not there is evidence of a correlation 
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between the presence of debris disks and close-in planets. 



4. Are Debris Disks and Close-in Planets Related Phenomena? 

Using the FEPS data, we address the possibility of a debris/planet connection by com- 
paring the results for the nine planet-bearing stars (hereafter "planet sample") to those 
of a larger sub-set of stars in the FEPS sample without regard to the presence of planets 
(hereafter "control sample"). The planet sample is a subset of the control sample. Given 
the current statistics from RV surveys, it is unlikely that the majority of the stars in the 
control sample harbor a giant planet, therefore, the control sample is likely less biased for 
the presence of giant planets than is the planet sample. By comparing these two samples, we 
investigate if the frequency and luminosity of debris disks are correlated with the presence 
of a massive planet. 



4.1. Selection of the control sample 

The main criterion for the choice of the control sample is that the observations reached 
similar levels of sensitivity as for the planet sample. For this we require that: (1) the stars 
in the control sample span the same range of distances as the planet bearing stars (26-62 
pc); and (2) their infrared background levels at 70 fim are similar. Age may also be a factor, 
as the stars in the FEPS sample range from 3 Myr to 3 Gyr typically, with a few stars 
perhaps as old as 3-10 Gyr, while planet bearing stars are typically older than 1 Gyr. If 
debris disks evolve significantly over the range 3 Myr-10 Gyr, this could also introduce a 
bias in the debris disk detection. Samp les of young stars show an initial rapid decline over 
the first ~100 Myr, while Bryden et al. (120061 ) found that the excess rate at 70 fim is 13±5% 



(down to a fractional luminosity of L^ ust /L*~10~ 5 , i.e. about 100 times the luminosity of 
the Kuiper Belt dust) and is not correlated with stellar age on Gyr timescales. Therefore, 
by choosing a control sample that is restricted to stars older than 300 Myr, we do not 
expect to introduce any significant age bias, while improving the statistics by increasing 
the number of stars in the control sample. Our control sample thus consists of 99 stars 
with distances 26-62 pc an d ages >300 Myr. We use the Kolmogorov-Smirnov (K-S) test 



(see e.g. Press et al. Il993l ) to assess whether or not the distributions of distances and IR 
background levels of the planet and the control samples are consistent with having been 
drawn from the same parent population. The K-S test yields two values: D, a measure 
of the largest difference between the two cumulative distributions under consideration and 
Probability(D > observed), an estimate of the significance level of the observed value of 
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D as a disproof of the null hypothesis that the distributions come from the same parent 
population, i.e. a very small value of Probability(D > observed) implies that the distributions 
are significantly different. Because in this case we find ProbabilityiD > observed) = 0.6 (for 
distance) and Probability (D > observed) = 0.4 (for IR background), we conclude that both 
samples could have been drawn from the same distribution in terms of distance and IR 
background levels and therefore can be compared. Note that the K-S test for age yield a 
much lower probability (Probability (D > observed) ~10~ 5 ), i.e. both samples are likely not 
drawn from the same distribution in terms of age. However, given that the observations 
indicate that for the ages under consideration (> 300 Myr, with approximately half of the 
stars having ages > 1 Gyr) there is no correlation between the 70 /im excess and the stellar 
age, we do not expect to introduce any significant bias by comparing both samples (but keep 
in mind that the validity of the comparison relies on the observed lack of correlation with 
age). 



4.2. Frequency of Debris Disks 

With respect to the frequency of debris disks, we find that 1/9 stars in the FEPS planet 
sample has 70 /im excess emission with a SNR in the ex cess >3, compared to 9/99 stars in 



the FEPS control sample; for the Bryden et al. (120061 ) survey, the rates are 1/11 (planet 
sample) and 7/69 (control sample). Because the frequency of debris disks (seen at 70 //m) in 
the planet sample and the control sample are similar, we conclude that there is no evidence 
of the presence of a correlation between the frequency of debris disks and close-in planets (if 
we were to assume a vA error in the number of stars with excesses, the frequency of debris 
around a planet-bearing star would be within a factor of 3 of the control sample). At 24 /im 
(tracing warmer dust), the frequency of debris disks could also be similar in the planet and 
the control samples, as none of the stars in the FEPS planet sample show excess emission, 
while 2/99 stars in the FEPS control sample do. 



4.3. Fractional Excess Luminosity: Survival Analysis 



The planet sample and the control sample are dominated by upper limits, therefore, the 
K-S test is not sufficient to assess the probability that they could have been drawn from the 
same parent distribution. To extract the maximum amount of information from the non- 
detections it is necessary to use survival analysis methods, which make certai n assu mptions 
about the underlying distributions. Using ASURV Rev 1.2 (LaValleyet al. Il992l ). which 



implements the survival analysis methods of Feigelson & Nelson (119851 ). we have used the 
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Gehan, logrank and Peto-Prentice tests to compute the probability that the planet sample 
and the control sample could have been drawn from the same parent distribution with respect 
to the fractional excess luminosity, L^t/L^. 

We use the fractional luminosity of the excess, Ld U st/L±, instead of t he 70 /im excess 



flux to minimize any correlation with distance. Following Bryden et al. (120061 ). from the 
70 /im excess emission one can estimate the fractional luminosity of the excess by assuming 
a single dust temperature, T dltst = 52.7 K, corresponding to an emission peak at 70 /im. 
In this case, L dust /L* ~ 10~ 5 (5600/T*) 3 (F 7CMust /F 70;Jl ,), where F 704ust and F 70i * are the dust 
excess and photospheric flux at 70 /im and T* is the stellar temperature in units of Kelvin 
. For non-detections, F 70t d us t = 3x AF 70 , where AF 70 is the 1-er uncertainty of the observed 
flux. 

The resulting survival analysis probabilities, using 3-cr upper limits, are 0.6 4 (Ge han), 



0.86 (logrank) and 0.72 (Peto-Prentice). As discussed in Feigelson & Nelson (119851 ). the 
logrank test is more sensitive to differences at low values of the variable under consideration 
(i.e. near the upper limits), while the Gehan test is more sensitive to differences at the 
high end (i.e. for the detections). The Peto-Prentice test is preferred when the upper limits 
dominate and the sizes of the samples to be compared differ (as it is our case). Similarly, we 
hav e carr ied out survival analysis for the sample of 69 FGK main sequence stars in Bryden et 



al. (120061 ). This sample was selected with regard to expected signal-to-noise ratio for stellar 
photospheres and is not biased for or against known planet-bearing stars. The planet sample 
consists of 11 stars with known close-in planets and the control sample includes all 69 stars. 
With respect to the FEPS targets, these stars are generally closer and the observations are 
therefore sensitive to less luminous debris disks (see Fig. 3). In this case, the probabilities 
that the planet sample and the control sample could have been drawn from the same parent 
distribution with respect to the fractional excess luminosity are 0.83 (Gehan), 0.86 (logrank) 
and 0.70 (Peto-Prentice). If we consider the FEPS and Bryden's samples together, these 
probabilities are 0.62, 0.85 and 0.70, respectively. Because all the probabilities are larger 
than 0.6, i.e. significantly larger than 0, the conclusion from the Gehan, logrank and Peto- 
Pren tice te sts from the data collected so far (from both FEPS and the GTO results in Bryden 



et al. 120061 ). is that we cannot rule out the hypothesis that the planet sample and the control 
sample have been drawn from the same population with respect to the fractional excess 
luminosity. In other words, we find no sign of correlation between the excess luminosity and 
the presence of close-in massive planets. 
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Discussion 



5.1. Comparison to Previous Studies 



Greaves et al. (12004 ) searched for submillimeter dust emission around 8 stars known 
from radial velocity studies to have giant planets orbiting within a few AU, and found no de- 
bris disks down to a dust mass limit of 6 x 10- 8 Mj; they also noted that out of 20 solar-type 
stars known to have disks, only one, e Eridani, has a planet orbiting inside a few AU (Hatzes 
et al. |2000| ). concluding that either debris disks and close-in giant planets are unrelated 
phenomena or they are mutually exclusive. However, these results had severe limitations 
due to the low sensitivity of the submillimeter observations and the limited accuracy of the 
radial velocity measurements due to the youth of the sample. The high sensitivity of the 
mid-IR observations with Spitzer allow a fresh reexamination o f the question. Preliminary 
work from a <Spzteer/MIPS GTO program by Beichman et al. ( 120051 ) found that out of 26 
FGK field stars known to have planets through radial velocity studies, six stars (HD 33636, 
HD 50554, HD 52265, HD 82943, HD 117176 and HD 128311) show 70 fjm excess with a 
SNR in the excess of 12.4, 12.1, 3.2 15.7, 3.2 and 4.6, respectively, implying the presence 
of cool dust (<100 K) located mainly beyond 10 AU. These six stars have a median age of 
4 Gyr and their fractional luminosities, l^duxt. I Lsfr y , in the range (0.1-1.2)xl0 -4 are ~100 
times that inferred for the KB (Beichman et al. 120051 ) . The study suggested a tentative 
correlation between the frequency and magnitude of the dust emissio n wit h the presence of 
known planets. Our analysis of the FEPS and the Bryden et al. (120061 ) samples do not 
confirm the presence of such a correlation. 



5.2. Interpretation of our Results 

We found that, given the Spitzer and radial velocity data we have so far, there is no 
evidence of a correlation between the presence of close-in massive planets and the frequency 
and excess luminosity of debris disks, i.e. debris disks are not more prevalent in systems 
with close-in massive planets than in systems selected without regard to whether they have 
close-in massive planets or not. 

This might seem a surprising result because it is reasonable to assume that most giant 



^Dus t mass estimates for the KB dust disk range from a total dust mass < 3x 10 10 M (Backman et 
al. ll995h to ~ 4xl0~ n M for dust parti cles < 150 (im (Moro-Martm & Malhotra liooj ) ; with a fractional 
luminosity of L ( / l(JS t/L*~10 -7 10 -6 (Stern 1996). The fractional luminosity of the asteroid belt dust (a.k.a 



zodiacal cloud) is estimated to be L,just/L*~10 8 -10 7 (Dermott et al. 2002h . 
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planets formed in systems that were initially rich in planetesimals, as plantesimals are the 
building blocks of giant planets in the core accretion model. However, despite a likely initial 
abundance of dust-producing planetesimals, systems with giant planets may not produce 
abundant debris at Gyr ages. The solar system is one example where there is significant 
evidence that it once had a massive planetesimal belt despite the little debris produced 
today. That is, giant planets may play an important role in the evolution of debris disks by 
efficiently grinding away or ejecting planetesimals from an initially massive disk. This could 
involve processes similar to the LHB in the early Solar System, where a large fraction of the 
dust producing planetesimals were lost due to the orbital migration of the giant planets. For 
the planet-bearing stars, given that the conditions for the formation of at least one giant 
planet were met, we could speculate that additional massive planets possibly formed and 
migrated, which could lead to LHB-type of events. Comparison of the debris disk properties 
between stars with and withou t mas sive giant planets may therefore be a function of age. 



The FEPS and Bryden et al. (120061 ) samples are drawn mainly from stars 300 Myr to 10 



Gyr old, i.e. mostly after the LHB is thought to have occurred in our Solar System. 

Our result also suggests that massive planets may not be required to produced debris. A 
possible mechanism for the production of debris in the pr esence or absence of massive planets 



is the collisional cascade model of Kenyon & Bromley (120051 ). Such a model can produce 
debris at Gyr ages, even in disks that are too low in solids to form a giant planet (i.e. too 
low in initial disk mass and /or too low in metallicity). In this model, large 1000 km size 
planetesimals can stir up smaller planetesimals (0.1-10 km in size) along their orbits, starting 
a collisional cascade that can produce dust excess emission of the magnitude shown in Fig. 
3 over the relevant range of ages. However, this cannot be the only mechanism because 
if it were to dominate debris production one would expect to see the dust tempera ture t o 



be correlated with age and this trend has not been observed (Najita & Williams 120051 ). 
Similarly, the observations in Fig. 3 could not confirm the time dependence of the fractional 
70 /im excess emission predicted by the models. 

That massive planets may not be required to produce debris is also supported by several 
observational results. Firstly, debris disks are more common than massive planets: it is found 
that >7% of stars have giant planets with M<13 Mj up and semimajor axis within 5 AU, 
but this is a lower limit because the duration of the surveys (6-8 years) limits the ability 
to detect planets between 3 AU and 5 AU. The expected frequency of gas g iant p lanets 
increases to ~12% when RV surveys are extrapolated to 20 AU (Marcy et al. I2005L with 
the distribution of planets following dN/dM ocM~ im from Mg atm , n to 10 M Jup (the surveys 
are incomplete at smaller masses). In comparison, t he fre quency of debris disks observed 



at 70 (Ltm with Spitzer is 13±5% (from Bryden et al. 120061 ) . Howe ver, this detection rate is 



sensitivity limited because the observations in Bryden et al. ( 120061 ) can only reach fractional 
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luminosities of Ld ust /L i >10 5 , i.e. >100 times the luminosity from our Solar System Kuiper 



Belt. Bryden et al. (120061 ) found that the frequency of dust detection increases steeply as 



smaller fractional luminosities are considered, going from nearly 0% for L^t/L^lO -3 , to 
2±2% for L dust /L*~10~ 4 and 13±5% for L dMSt /L^~10~ 5 . Using this cumulative distribution 
and a ssuming that the distribution of debris disk luminosities is a Gaussian, Bryden et al. 
( 120061 ) estimated that the luminosity of the Solar System dust is consistent with being 10 
x brighter or fainter than an average solar-type star, i.e. debris disks at the Solar System 
level could be common. The debris disks observed with Spitzer could therefore be the high 
luminosity tail of a distribution of dust luminosities that peaks near the Solar System values. 

Secondly, there is no c orrela tion between stellar metallicities and the incidence of de - 
bris disks (Beichman et al. 12005k Bryden et al. 120061 and Greaves, Fischer & Wyatt 120061 ). 
Greaves, Fischer & Wyatt (120061 ) found that in a sample of 310 F7-K3 stars within 25 pc of 
the Sun and for which the stellar metallicities are known, there is only a 0.6% probability that 
planet-bearing stars and debris disks stars have the same metallicity distribution, w ith th e 
planet-bearing stars being correlated with high stellar metallicities (Fischer & Valenti 120051 ) . 
This is in agreement with the core accretion model, where the formation of giant planets 
requires the presence of a large surface density of solids in the disk, so that the planet can 
grow a core sufficiently large to accrete an atmosphere before the gas disk disappears in 
<10 Myr. Because the governing time scale in the growth of planetesimals is the orbital 
period, in the KB region the planetesimal formation process is slower (according to Kenyon 
and Bromley 12004 it may take ~ 3 Gyr to form a Pluto at 100 AU), but can proceed well 
after the gas disk has dissipated (so there is no time limitation). This can occur in systems 
regardless of whether or not they meet the conditions for giant planet formation. 

As a result, collisional grinding in a self-stirred model of this kind might be expected 
to produce debris in systems with low metallicities and low initial disk masses. If this leads 
to debris production in a wider variety of systems than can produce giant planets, we might 
expect the presence of debris to be poorly correlated with the presence of giant planets. 



6. Conclusions 

• Using Spitzer observations, we have searched for debris disks around 9 planet-bearing 
solar-type stars, with stellar ages ranging from 2 to 10 Gyr. Only one of the sources, 
HD 38529, is found to have excess emission above the stellar photosphere, with a signal- 
to- noise ratio at 70 fim of 4.7 and no excess at A < 30 /im. The remaining sources 
show no excesses at any Spitzer wavelengths. 

• Given the data we have so far, from both FEPS and the FGK sample from Bryden 
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et al. (120061 ). and using survival analysis, we find that there is no evidence of a 
correlation between the presence of close-in massive planets and the frequency and 
excess luminosity of debris disks. 

Because we expect massive planets to form in systems that are initially rich in plan- 
etesimals, but the observations indicate that systems with giant planets do not prefer- 
entially show debris, there is the possibility that massive planets play an important role 
in the evolution of debris disks by efficiently grinding away or ejecting planetesimals 
from an initially massive disk, possibly in a LHB-type of event. 

Our result also suggests that massive planets may not be required to produced debris , 



which is supported by the collisional cascade models of Kenyon & Bromley (120051 ). 
and the observations and theoretical models that indicate that debris disks are more 
prevalent than massive planets. 

APPENDIX: HD 150706 

HD 150706 is part of the FEPS sample and exhibits an excess emission at 70 /jm with a 
SNR in the excess of 4.3 (and a color-corrected flux of 46.3 mJy). Even though it has been 
listed as a planet-bearing star, it is not included in our planet sample because new radial 
velocity observations cannot confirm the claimed planet. HD 150706 has appeared in variou s 



compilations of sun-like stars with extra-solar planets (c.f. Santos, Israelian, Mayor 120041 ) . 
An orbital solution for a purported 1.0 Mj up eccentric planet at 0.8 AU was announced by the 
Geneva Extrasolar Planet Search Team (2002, Washingto n con ference "Scientific Frontiers in 



Research in Extrasolar Planets"; Udry, Mayor & Queloz. 120031 ); however, there is no refereed 
discovery paper giving details, only web pages. Eight doppler velocity measurements (Table 
4) made with HIRES on the Keck telescope from 2002 to 2006 yield RMS of 12.1 m/s, far 
below the 33 m/s claimed velocity amplitude due to a planet. The RMS for a linear fit of the 
HIRES data is 8 m/s which can be adequately explained by the expected jitter for a young 
(700±300 Myr) and active early G star like HD 150706. The four years of HIRES data rule 
out the presence of planets of roughly 1 Mj up or larger located within 2 AU, and 2 Mj up 
or more within 5 AU (modulo sin(i)). Smaller planets inward of 5 AU, and super- Jupiters 
outward of 5 AU are not inconsistent with the HIRES observations to date. Further, the 
lack of a monotonic trend in the velocities of amplitude many tens of m/s indicates that 
there is no brown dwarf nor a low mass star anywhere within ~ 20 AU. 
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Table 1. Stellar properties 



Source 


Spectral 


Distance a 


Age 


T e// 


log(L) 


M 


[Fe/H] 


(HD #) 


Type 


(PC) 


(Gyr) 


(K) 


(log(L )) 


(M ) 




6434 


G2/3V d 


40±1 


12±l c 


5835 s 


0.05±0.02 c 


0.84±0.05 f 


-0.52 s 


38529 


G8III/IV d 42±2 


3.5±l h 


5697 s 


o.so s 


1.47 s 


0.445 s 


80606 


G5J 


58±20 


6 k 


5573 s 


-0.15 s 


1.06 s 


0.343 s 


92788 


G6V 1 


32±1 


6±2 m 


5836 s 


o.or 


1.13 s 


0.318 s 


106252 


GOJ 


37±1 


5.5±l n 


5870 s 


0.11 s 


l.or 


-0.076 s 


121504 


G2V° 


44±2 


2±P 


6075 s 


0.19±0.04 c 


1.03±0.06 f 


0.16 s 


141937 


G2/3V q 


33±1 


2.6±l r 


5847 s 


0.07 s 


1.08 s 


0.129 s 


179949 


F8V q 


27±1 


2.5±l s 


6168 s 


0.27 s 


1.21 s 


0.137 s 


190228 


G5IV 


62±3 


5 s 


5348 s 


0.63' 


1.21 s 


-0.180 s 



Referee 
from FEPI 
et al. (200 
et al. (201 
Fischer (2i 
(1) Houk I 
Gonzalez < 
et al. (20( 
& Gratton 
Barbieri & 
Fischer (2i 
Fischer (2( 
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Table 2. Orbital characteristics of known planetary companions 



Planet 


M p sim 


Period 




6p 


Nobs 


Ref. 


(HD #) 


(M Jnp ) 


(days) 


(AU) 








6434b 


0.397(59) 


21.9980(90) 


0.1421(82) 


0.170(30) 


130 


(1) 


38529b 


0.852(74) 


14.3093(13) 


0.1313(76) 


0.248(23) 


162 


(2) 


38529c 


13.2(1.1) 


2165(14) 


3.72(22) 


0.3506(85) 


162 


(2) 


80606b 


4.31(35) 


111.4487(32) 


0.468(27) 


0.9349(23) 


46 


(2) 




3.90(9) 


111.81(23) 


0.47 


0.9227(12) 


61 


(3) 


92788b 


3.67(30) 


325.81(26) 


0.965(56) 


0.334(11) 


58 


(2) 




3.58 


325.0(5) 


0.96 


0.35(1) 


55 


(1) 


106252b 


7.10(65) 


1516(26) 


2.60(15) 


0.586(65) 


15 


(2) 




7.56 


1600(18) 


2.7 


0.471(28) 


40 


(4) 


121504b 


1.22(17) 


63.330(30) 


0.329(19) 


0.030(10) 


100 


(1) 


141937b 


9.8(1.4) 


653.2(1.2) 


1.525(88) 


0.410(10) 


81 


(5) 


179949b 


0.916(76) 


3.092514(32) 


0.0443(26) 


0.022(15) 


88 


(2) 


190228b 


4.49 


1146(16) 


2.25 


0.499(30) 


51 


(4) 



Table 3. Spitzer photometry and Kurucz stellar models for FEPS targets with planets 



Source 




IRAC 3.6 


IRAC 4.5 


IRAC 8 


IRS 13 


MIPS 24 


IRS 24 


IRS 33 


MIPS 70 


SNR e:rc 


1200 


(HD #) 




(/im) 


(H 


(H 


(/im) 


(/im) 


(/im) 


(/im) 


(/im) 


at 70 /xm 


(/im) 


d A O A 

b4o4 


obs 


n r o i 7 

yoz± f 


o0o± ( 


Zlo±l 


7/1 n_i_n o 
f4.y±0.o 


oo n_t_n o 
Zo.y±0.z 


o c n_i_n o 
ZD. 0±0.O 


1 1 _i_ 1 
11±1 


Q n_l_7 A 


n 7 
0. / 


n_i_i n 
0±10 




model 


oyZ±Zo 


000±1 ( 


z(Jo±o 


/ o±Z 


no n i n 7 


no o i n 7 
Zo.Z±U. / 


i o o_i_n /! 
1Z.Z±U.4 


Z.OztU.l 






no tron 

ooozy 


obs 


tronQ-L/) O 

ooyo±4z 


O^O /I _1_ A A 

obo4±44 


1 o /i n_i_n 
lo40±y 


A C 7_1_ K 

4o I ±o 


1 txn_i_ 1 

lo0±l 


1 /I £ 1 o 

14b±z 


oo±z 


7r i i -1 


/I 7 

4. / 






model 


Oyo0±zoo 


dO 

oooy±loo 


loo0±oo 


4o l ±Z4 


i k£?_i_q 
l00±O 


i r /? i o 

loo±o 


oz±4 


i 7 i _i_n o 
1 ( .4±0.o 






oOoOo 


obs 


o a n_i_o 
o40±Z 


1 n i o 
Z10±O 


/ o.o±0.o 


no Q_l_n 7 

Zo. o±0. / 


o r? r in no 

o.oo±0.0o 


7 n_t_n o 


o 1 _i_n £ 
O.l±0.O 


o i r 
C)±0 


n c: 
0.0 






model 


66o±lZ 


OP1SX7 

ZUo± / 




Z f ±1 


o.o±U.o 


o.o±U.o 


4.0±U.Z 


U.yo±U.U4 






yz ( so 


obs 


144 / ±1U 


oyl±l 1 


6Z6±Z 


11111 

11 1±1 


Q£? 1 _i_n q 


97 11 


on i o 
ZU±Z 


1 1 i n 

1 l±y 


U.o 


r i 1 r 
D±lD 




model 


1 /I io i ro 

14oo±0Z 


nno_i_o 1 
y0z±ol 


oon_i_ i o 
ozy±lZ 


1 1 o_i_ /I 

llo±4 


QO I 1 

oo±l 


OO I 1 

oo±l 


1 n o _i_n 7 
ly.O±0. I 


a o_i_n o 
4.Z±0.2 


1 




lUbZoZ 


obs 


lZUU±y 


/ 4o±y 


Z ( 1±Z 


no o_i_n n 
yz.o±U.y 


on p.-L.c\ o 
oU.0±U.o 


oU.O±U.O 


~\ O I 1 

lo±l 


lo±y 


1 /l £3 
1.4 o 






model 


T T /1£_i_Q1 

114o±ol 


/ Uo±lo 


Q/?Q 1 7 

Zoo± / 


O/l 1 o 

y4±o 


on n_i_n q 
o(J.U±U.o 


on n_i_n q 
oU.U±U.o 


l0.o±U.4 


o.o±U.l 


1 




121504 


obs 


1002±7 


631±8 


225±1 


81.3±0.8 


25±0.2 


27±2 


14±2 


27±19 


1.3 






model 


976±35 


620±22 


222±8 


79±3 


25±1 


25±1 


13.4±0.5 


2.8±0.1 






141937 


obs 


1393±10 


872±11 


311±2 


110±1 


34.9±0.3 


35.1±0.5 


18±2 


-3±12 


-0.6 






model 


1365±38 


864±23 


312±9 


111±3 


36±1 


36±1 


18.7±0.5 


4.0±0.1 






179949 


obs 


2943±21 


1849±23 


658±4 


234±2 


73.9±0.7 


71±1 


35±2 


-5±11 


-1.2 


1±10 




model 


2822±78 


1809±49 


641±18 


229±6 


73±2 


73±2 


38±1 


8.1±0.2 






190228 


obs 


2068±15 


1283±16 


469±3 


166±2 


52.8±0.5 


53.4±0.8 


29±3 


12±26 


0.2 






model 


1987±66 


1223±40 


456±15 


164±5 


52±2 


52±2 


27.6±0.9 


5.8±0.2 







NOTE. - Photometry and 1-a internal uncertainties are in units of mJy. Calibration uncertainties are not included in the error 
estimates. IRS fluxes come from synthetic photometry from IRS low-resolution spectra. Fluxes at 1200 /im are from Carpenter et al. 
(2005). (model) is the expected stellar contribution from its Kurucz model; (obs) is the photometric measurement. The signal-to-noise 
ratio of the excess, SNR e:rc , is the photometric measurement minus the star's contribution from its Kurucz model divided by the global 
uncertainty. The global uncertainty is calculated adding in quadrature the internal and calibration uncertainties, the later taken to 
be 7% for MIPS 70 pim. 
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Table 4. Relative Radial Velocities for HD 150706 



JD 


RV 


Unc. 


-2450000 


(m/s) 


(m/s) 


2514.728 


-0.40 


1.8 


2538.717 


-3.59 


1.9 


2713.091 


4.20 


1.7 


2806.008 


-9.49 


2.0 


2850.896 


-5.08 


1.8 


3179.924 


0.00 


1.7 


3427.128 


23.83 


2.2 


3842.061 


20.89 


2.2 
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Fig. 1. — Spectral energy distributions (SEDs) of the nine planet-bearing stars in the FEPS 
program. The green line is the Kurucz model. The black thicker line is the IRS low-resolution 
spectrum. The photometric points are identified as follows: green squares are ground based 
observations (including Tycho and 2MASS); black circles are Spitzer observations (IRAC, 
MIPS and synthetic photometry from IRS); red diamonds are IRAS observations. In all cases, 
the error bars correspond to l-a uncertainties. Upper limits are represented by triangles and 
are given when F/AF < 3 and placed at F + 3xAF if F > 0, or 3xAF if F < 0. Black 
triangles are upper limits for Spitzer 70 urn, and red triangles for IRAS and 1.2 mm. 
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Fig. 2. — Color-color diagrams of the 9 stars in the FEPS planet sample (left) and the 99 
stars in the FEPS control sample (right). Stars that show 70 /im excess emission with a 
SNR in the excess > 3 are shown in red, and include one star in the planet sample (HD 
38529) and nine stars in the control sample. Similarly, stars that show 24 /im and 33 jum 
excess emission with a SNR in the excess > 3 are shown in blue and green, respectively, and 
include only two stars (for 24 /im) and four stars (for 33 /im) in the control sample and none 
in the planet sample. [The outlier at (0.114, 0.357) in the upper left panel corresponds to 
HD 80606: its IRS spectrum is very noisy beyond 30 microns, possibly due to a peak-up 
failure, making the 33 /im point unreliable]. 
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Fig. 3. — Ratio of the excess flux to the photospheric flux for the stars with 70 /im excess 
emission and a SNR in the excess > 3. The shape of the symbol indicates the presence 
of a close-in planet. Circles are stars with known radial velocity planets; crosses are stars 
without known planets (but included in radial velocity surveys); and triangles: stars without 
known planets (not included in radial velocity surveys). The black symbols corresp ond to 
stars in the FEPS survey , while blue and red correspond to stars in Beichman et al. (120051 ) 
and Bryden et al. (120061 ). respectively. (The latter includes stars at smaller distances than 
those in the FEPS sample, so is sensitive to smaller excesses). For comparison, the dotted 
line shows Fd ust /F* at 60 /im res ulting from the collisional cascade of a planetesimal disk 
at 30-80 AU (Kenyon & Bromley 120051 ). Because the collisional physics and the behavior 
of the debris following the collision are uncertain, the estimate for Fd usf /F* at 60 /im could 
vary by more than a factor of 10, so the observations could be consistent with the model 
predictions. 



